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ABSTRACT 

The standard structure formation model based on a ACDM cosmology pre- 
dicts that the galaxy clusters have triaxial shapes and that the cluster galaxies 
have a strong tendency to be located preferentially along the major axes of host 
cluster's dark matter distributions due to the gravitational tidal effect. The 
predicted correlations between dark matter and galaxy distributions in triaxial 
clusters are insensitive to the initial cosmological parameters and to the galaxy 
bias, and thus can provide a unique test-bed for the nonlinear structure forma- 
tion of the ACDM cosmology. Recently, Oguri et al. determined robustly the 
dark matter distributions in the galaxy clusters using the two dimensional weak 
lensing shear fitting and showed that the orientations of the cluster galaxy distri- 
butions are only very weakly correlated with those of the underlying dark matter 
distributions determined robustly, which is in contrast to with the ACDM-based 
prediction. We reanalyze and compare quantitatively the observational result 
with the ACDM-based prediction from the Millennium Run simulation with the 
help of the bootstrap resampling and generalized x 2 -statistics. The hypothesis 
that the observational result is consistent with the ACDM-based prediction is 
ruled out at the 99% confidence level. A local fifth force induced by a non- 
minimal coupling between dark energy and dark matter might be responsible 
for the observed misalignments between dark matter and galaxy distributions in 
triaxial clusters. 

Subject headings: cosmology:theory — large-scale structure of universe 



1. INTRODUCTION 



Although the dark energy that is responsible for the present cosmi c acceleration is know n 
to occupy more than 70% of the energy budget of the Universe (e.g., Komatsu et al.ll2010l ). 
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its nature is still shrouded in deep mystery. The simplest candidate for the dark energy is 
the cosmological constant (A) that Einstein int roduced in his field equation a century ago 
and called the fundamental constant of Nature ( IEinsteinlll917l ). which has been interpreted 
as the vacuum energy that is completely homogeneous in space, does not evol ve with time 



non-i nteracting;, chemically inert, and has a constant equation of state (see ICarroll et al. 
1992, for a review). 



Over the past two decades, the cosmologists have witnessed and enjoyed the practical 
success of the ACDM (cosmological constant + cold dark matter) cosmology in explaining 
the large-scale properties of the Universe, which have left little doubt on the validity of the 
ACDM as a standard model. Nevertheless, an increasing number of literatures have been 
devoted to seeking for viable alternatives that could compete with and hopefully replace the 
ACDM (lAlcanizI 12006k for a review). The main motivation of the alternative dark energy 



models is to explain naturally the cosmic coincidence that the densities of dark energy and 
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An intriguing idea that was motivated by a fundamental particle physics and has gained 
a lot of popularity recently in the field of large scale stru cture is that dark energy is a scalar 
field having a non- minimal coupling with dark matter (lAmendolal l2000t iFarrar fc Peebles 
2004 ; lAlcanizI 120061 and references therein). The distinct feature of this fascinating idea is 
that a possible fifth force can be generated by the propagation of dark matter inhomogeneities 
into the scalar field due to the non-minimal coupling. The strength, range and potential 
form of a possible fifth force depends on how to model the coupling in the dark sector (dark 
energy+dark matter). The reason for the recent sharp attention on the coupled scalar field 
models is that the inclusion of a fifth force generat ed by the dark sec t or interaction may 



be able to explain seve r al observational p henomena ( jMaccio et al.l 12004 ; iNusser et al.l 12005 
Peebles fc Nusserll2010l ; iBaldi et al.ll2010l ). which have been suggested as mismatches to the 
theoretical predictions of the ACDM cosmology such as th e abundance of ga l axy satellites 



the speed of bullet cluster, and the clear-out of void dwarfs ( IKlypin et al.lll999t iPeeblesI 12001 
Farrar fc Rosen] 120071 : iLee fc Komatsulboioh . 



In this Letter, we present another observable that exhibits an obvious mismatch with 
the numerical prediction based on the standard structure formation based on a ACDM cos- 
mology: the misalignments between matter and galaxy distributions in triaxial galaxy clus- 
ters. In conventional approaches the shapes of dark matter distribution in galaxy clusters 



were often determined by measuring the member galaxy distributions (e.g.. iBinggeli 
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Evans k, Bridle! 2009 ). Recent progress in weak lensing shear 
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analysis, however, has allowed us to determine the shap es of clusters directly from dark mat- 
ter distributions (jOkabe et al.ll2010l ; lOguri et al.ll2010l ). which has revealed that the cluster 
dark matter distributions are only very weakly aligned with the cluster galaxy distributions. 

Noting that the standard structure formation model based on ACDM cosmology pre- 
dicts very stron g correlations between the shapes of galaxy and dark matter distributions 
( ILee et al.l 120081 ). we make a quantitative comparison between the observed trend and the 
ACDM-based prediction, and make a robust test of the hypothesis that the observed trend 
is consistent with the ACDM-based prediction with the help of the bootstrap resampling 
and generalized x 2 statistics. We speculate that the observed misalignments between dark 
matter and galaxy distributions in triaxial clusters might be a new observational signature 
of a dark sector interaction. 



NUMERICAL RESULT BASED A ACDM COSMOLOGY 



The correlations between the major ax es of dark matter and satellite galaxy distributions 
in dark halos have been investigated by iLee et al.l (120081 ) using the data from the mili- 
Millennium simulation that was run on a periodic box size of 62.5 h~ l Mpc fo r a flat ACDM 
cosm ology with Q m = 0.25, = 0.75, h = 0.73, cr 8 = 0.9, n s — 1 (jSpringel et al. 
20051 ). The dark halos and their substructures were identified from the mili-Millennium 
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semi- analytic galaxy catalog, we refer the readers to ISpringel et al.l (120051 ) . 



Selecting those mili-Millennium halos which have five or more satellite galaxies, ILee et al. 



(120081 ) ended up having a sample of 277 dark halos with mass M > 1.3 x 10 13 h~ l M Q at 
z = 0. They determined the major axes of the selected 277 halos using two different schemes. 
The first scheme is based on the dark matter distribution of each halo. They calculated the 
inertia momentum tensor of each halo using the constituent dark matter particles (each of 
which has mass 8.6 x 1O 8 /i _1 M ) and determine the major axis as the eigenvector of the 
inertia momentum tensor corresponding to the largest eigenvalue. In the second scheme the 
major axis of the inertia momentum tensor of each halo was calculated by using the satellite 
galaxies weighted by their luminosity. 



The correlations between the major axes determined by the two different schemes were 
quantified in terms of the distribution of the alignment angle: Basically, they measured the 
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angle between the two major axes in range of [0, 90°] for each halo. Binning the alignment 
angles and counting the number of halos whose alignment angles belong to each bin, they 
calculate the number fraction distribution as a function of the angle between the two major 
axes of dark matter and satellite galaxy distributions (see Figure 1 in Lee et al. 2008). Very 
strong correlations between the two major axes were found, which proved that the satellite 
galaxies in the massive halos tend to be located near the major axes of the dark matter 
distributions. 

The correlations between matter and galaxy distributions in halos are usually ascribed 
to the effect of the exter nal tidal fields: A ccording to the cosmic web theory based on a cold 
dark matter paradigm (IBond et al.lll996l ). the coherence and non- linear sharpening of the 
gravitational tidal fields induce the shape-alignments between the large- scale structures such 
as cluster-supercluster alignments, v oid-supercluster alignments and the inclination of the 
galaxy spin axes onto the sheets (e.g., 



Navarro et al. 



20061 : iTruiillo et all 12009 : iLee & Evrardl 120071 Park fc Led 120071 1 . 



2004 



Kas un fc Evrardll2005l : I Alt ay et al. 



The gravitational tidal effects also cause the galaxy clusters to have triaxial shapes 
and induce the correlations between matter and satellite galaxy distribution in the triaxial 
clusters since the merging and accretion of matter and galaxies onto the clusters occur 
preferentially along the primary filaments that are aligned with the directions of t he minimal 



matt e r compression, i.e., the minor principal axes of the gravitational tidal fields (I West et al. 
19951 : iBailin et al.l 120051 ; iBailin fc Steinmetzl 120051 ) . After anisotropic merging along the 



primary filaments, however, the secondary infall and subsequent nonlinea r evolution insid e 
the clusters would modify the sp atial distributions of the member galaxies iDubinskil (119991 ). 
For instance, ILee fc Kangi (120061 ) have shown that as the clusters become more relaxed, the 
correlations between the principal axes of dark matter and substructures tend to decrease. 
They also found that the correlations increase with redshift and cluster mass. 

In observations, the correlations between the major axes of dark matter and galaxy 
distributions are often measured in the two dimensional projected planes. Therefore, for a 
direct comparison with the observational result, it is necessary to predict numerically the 
correlations of the projected major axes of dark matte r and sate l lite g alaxy distributions. 
Adopting the same mili-Millennium Run data used by ILee et al.l ( 120081 ) and assuming the 
flat-sky approximation, we investigate here the correlations between the projected major 
axes of dark matter and galaxy distributions. Measuring and binning the angles <j) of the 
projected major axes of dark matter and satellite galaxy distributions of the selected 277 
halos in the two dimensional plane, we find the number fraction of dark halos as a function 
of 6. 



The left, middle and right panel of Figure [T] plots the results for the case that the 
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projection was done onto the xy, yz, and zx plane, respectively. The errors include both of 
the sample variance as well as the Poisson noise. The sample variance is calculated as one 
standard deviation among the 1000 bootstrap resamples. As one can see, the distributions 
has a sharp peak at the first bin 0° < <fi < 10°, as in the three dimensional case, which 
proves that projections onto the two dimensional planes do not diminish the correlations 
between the major axes of matter and satellite galaxy distributions. In the following section, 
we compare this numerical prediction with the observational result and test the hypothesis 
that the observational result is consistent with the numerical prediction based on a ACDM 
cosmology. 



2.1. NUMERICAL VS. OBSERVATIONAL 



Performing two dimensional a nalysis of the weak le n sing shear maps co nstructed with 
high-quality Subaru/Supreme-cam (IMiyazaki et al.ll2002l ). IOguri et al.l (120101 ) have measured 
the two dimensional projected shapes and orientations of dar k matter d istributions in the 



galaxy clusters from the Local Cluster Substructure Survey (jOkabe et al.l 120101 ) . Using a 
subsample of 18 galaxy cl usters for which th e weak lensing shear maps are well fitted to the 
elliptical density profiles ( iJing fc Sutol 120021 ) . they determined the longest axis of the dark 
matter distribution and measured its position angle for each cluster (see Table 1 in Oguri et 
al. 2010). 

Selecting those member galaxies whose r-band magnitudes are brighter than 22 mag, 
they also fitted the member galaxy distribution (confined to the same region) in each cluster 
to an elliptical power law profile and determined the positions angle o f the longest axis o f 
galaxy distributions (see Table 2 in Oguri et al. 2010). According to lOguri et al.l ( 120101 ) . 
the fitting of the dark matter and member galaxy distribution in each cluster was done at 
the square region 20' x 20', around cluster's center chosen to be the location of the brightest 
cluster galaxy, where 20' corresponds to approximately 2-3 Mpc for the selected clusters in 
a redshift range of 0.1 < z < 0.3. For a detailed desc ription of th e two dimensional fitting 
of weak lensing shear field and the cluster sample, see lOguri et al.l ( 120101 ). 



By comparing the position angles of dark matter and galaxy distributions in each cluster, 
they found only very weak correlations between the two. They considered several possible 
systematics such as galaxy bias, dilution effect, fitting area and field-galaxy contamina- 
tion but found no systematics significant enough to explain the observed trend. Although 
Oguri et al.l ( 120101 ) asserted that the ellipticity distribution of dark matter distributions in 
the observe d galaxy clusters a grees well with the theoretical prediction based on a ACDM 
cosmology (IJing fc Sutol 120021 ) . they did not recognize that the observed misalignments be- 
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tween the longest axes of dark matter and member galaxy distributions in the clusters are 
inconsistent the ACDM prediction. 

To make a parallel comparison of the observational result with the ACDM prediction 
given in §2, here we determine the distributions of the angles between the projected major 
axes of dark matter a nd ga laxy distributions using the observational data of the 18 clusters 



given m 



Oguri et al.l ( 120101 ). The angle 4> for each cluster is determined by subtracting the 
position angle of the galaxy distribution from that of dark matter distribution and taking 
its absolute value. If the position angle difference exceeds tt/2, then we take n — (f>. Then, 
we bin the values of <fi and count the number of the clusters belonging to each bin. 

Figure [2] plots the result as open squares. The dotted line represents the uniform distri- 
bution corresponding to the case that there is no correlation between the major axes of dark 
matter and galaxy distributions. The errors include the sample variance and the Poisson 
noise. As in §2, the sample variance is calculated as la scatter among the 1000 bootstrap 
resamples. For comparison, we also plot the numerical prediction of the ACDM cosmology 
obtained in §2 as thick solid line (projection onto xy plane) while the thin solid lines represent 
the ±1(7 bootstrap errors of the numerical results. 

As it can be seen, the observational result deviates significantly from the numerical result 
based on a ACDM prediction and the semi-analytic galaxy formation model. To quantify the 
deviation, we employ the generalized x 2 -statistics which accounts for the correlations between 
the neighbor 0-bins due to the errors in the measurement of cluster's position angles: 

x 2 = >:,A//VA./^ (i) 

Here is the covariance matrix defined as Cy = (AfiAfj) and Afi is the difference 
between the observed and predicted value of the cluster number fraction at the i-th bin 
given as Afi = f° bs — f the (4>i). The bracket in the definition of (Equation [T] represents 
the ensemble average over the 1000 bootstrap resamples. Noting that in the large (f) section 
(beyond the third bins) the Poisson noises dominate over the bootstrap errors and the signals 
of the numerical results are confined to the first two bins, we consider the first three bins 
(0° < <p < 30°) for the calculation of x 2 (i- e -> the degree of freedom is 3). 

Putting the numerical results obtained in §2 into f the {<pi) for the calculation of x 2 , we 
test the hypothesis that the observational result is consistent with the ACDM prediction and 
found that the hypothesis is rejected at the 99% confidence level. Meanwhile, putting the 
uniform distribution into f the (4>i), we also test the null hypothesis that there is no correlation 
between the major axes of dark matter and galaxy distributions and have found that the 
null hypothesis is rejected only at the 23% confidence level. 



Excluding those clusters whose position angles suffer from large uncertainties, we end up 
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with a smaller subsample of 13 clusters and redetermined the number fraction distribution 
as a function of <fi. Figure [3] plots the results. As it can be seen, the distribution has a 
higher peak at the second bin but still deviates significantly from the numerical result based 
on a ACDM cosmology and the semi-analytic galaxy formation model. The recalculation of 
X 2 for this case rejects the hypothesis that the observational result from the 13 clusters is 
consistent with the numerical result is rejected at the 95% confidence level. 

It is worth mentioning the differences in the redshift range between the numerical and 
observational data used for the above comparison: The numerical result has been obtained 
at z = while the observational results have been drawn from the galaxy clusters at 0.1 < 
z < 0.3. This difference in the redshift range, however, would even worsen the disagreement 
between the numerical and the observational re sult for the following reason. As mentioned 



in §2 and as shown by N-body simu lations (e.g.. iBailin &: Steinmetzll2005l ; iBailin et al.ll2005 
Altay et al.ll2006t iLee fc Kanp2006l ). the correlations between the major axes of dark matter 
and satellite galaxy distributions tend to be stronger at higher redshifts. In other words, the 
predicted strength of the correlations used for the comparison with the observational result 
is underestimated. 

One might think that the selection bias in the measurements of galaxy distributions 
should be responsible f or th e disagreement between theory and observation since in the 
analysis of lOguri et al.l ( 120101 ) only those bright galaxies with mag > 22 at r-band are used 
unlike in the numerical analysis. Recall, however, that we used the luminosity-weighted 
galaxy distributions in the numerical analysis to determine the major axes, which should 
minimize the expected selection bias. 

Another difference lies in the values of the cosmological parameters used for the mili- 
Millennium Run simulations. Especially, the value of erg chosen b y the Millennium Ru n 
simulation has been known to be higher than the WMAP value (IKomatsu et al.l 120100 . 
The correlations between dark matter and satellite galaxy distributions in triaxial clusters, 
however, are unlike to be significantly affected by the initial cosmological conditions, since 
it represents a nonlinear observable rather than a linear one. 



3. DISCUSSION 

Now that we have found the observed misalignments of the projected major axes of dark 
matter and galaxy distributions in the galaxy clusters to be inconsistent with the numerical 
prediction based on a ACDM cosmology and the semi-analytic galaxy formation model, we 
would like to discuss the possible physical origin of the misalignments. A usual suspect would 
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be some hydro dynamic gas physics involved in the galaxy formation that was not included 
in the Millennium Run simulation and semi-analytic galaxy catalog. If this is really the case, 
then our result indicates that the semi-analytic galaxy formation model has to be improved 
to account for the observed deviations of the orientations of cluster galaxy distributions 
relative to that of dark matter distributions. 

An intriguing new possibility is that the orientations of member galaxy distributions 
have deviated from that of dark matter distributions due to the existence of the fifth force 
that is non-zero only in massive halos. In some dark energy model where the scalar field 
is coupled both to the dark matter particles and the second relativistic particle species, a 
fifth force is short ranged, h aying non-zero values only on the scales sm aller than the typical 



cluster sizes (~ IMpc) (e.g., Farrar &: Peebles! 2004 ; Nusser et al. 2005 ) 



According t o recent simulation results (e.g., iNusser et al.l 120051 ; iKeselman et al.l 12009 
Baldi et al.ll2010l ). in these models with a local fifth force the cluster galaxies (i.e., the galac- 
tic halos in dense environments) tend to form earlier and evolve more nonlinearly than in 
the standard ACDM model. In addition, a local fifth force also tend to affect the mutual 
interactions between the satellite galaxies, which would result in an internal tidal field dif- 
ferent from the gravitational external tidal field. Whereas, the dark matter particles which 
constantly accrete onto the host clusters at relatively later epochs than the member galaxies 
along the primary filaments from the surroundings would be less affected by the existence 
of a local fifth force and thus would keep the memory of the external gravitational tidal 
fields. Overall, a fifth force would play a role of diminishing the expected strong alignments 
between dark matter and galaxy distributions in clusters. 

This new idea also suggests that it might be possible to constrain the range and strength 
of a local fifth force by measuring the misalignments between dark matter and galaxy dis- 
tributions in clusters. This goal calls for making a quantitative prediction of the coupled 
dark energy model on the correlations between dark matter and galaxy distribution using 
detailed N-body and hydrodynamic simulations. 



The Millennium Run simulation used in this paper was carried out by the Virgo Super- 
computing Consortium at the Computing Centre of the Max-Planck Society in Garching, 



which are available at |http: / / www.mpa-garching.mpg.de / millennium We thank M.Oguri 



for useful comments. This work was supported by the National Research Foundation of 
Korea (NRF) grant funded by the Korea government (MEST, No.2010-0007819). 
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Fig. 1. — The fraction of the dark halos as a function of the angles between the major axes 
of the cluster halos and the spatial distributions of the member galaxies projected into the 
x-y, y-z plane and z-x plane (in the left, middle and right panel, respectively). 
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Fig. 2. — Comparison of the A prediction (solid line) with the observational result (open 
squares)obtained from Oguri et al. (2010). The error-bars include the sample variance as 
well as the Poisson noise. The thin solid lines represent the ±la bootstrap scatters. 
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Fig. 3. — Same as Figure |2] but including in the observational results only those clusters 
which do not show larger uncertainties in the measurement of the major axis directions. 



